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^<\ Abstract 

f~^ We study the scattering in a quantum star graph with a Fulop-Tsutsui coupling in its vertex 

^SJ and with external potentials on the lines. We find certain special couplings for which the prob- 

5_i ability of the transmission between two given lines of the graph is strongly influenced by the 

^ potential applied on another line. On the basis of this phenomenon we design a tunable quan- 

^> tum band-pass spectral filter. The transmission from the input to the output line is governed by 

a potential added on the controlling line. The strength of the potential directly determines the 
passband position, which allows to control the filter in a macroscopic manner. Generalization of 
this concept to quantum devices with multiple controlling lines proves possible. It enables the 
construction of spectral filters with more controllable parameters or with more operation modes. 
In particular, we design a band-pass filter with independently adjustable multiple passbands. We 
also address the problem of the physical realization of Fulop-Tsutsui couplings and demonstrate 
'S that the couplings needed for the construction of the proposed quantum devices can be approxi- 

C^ mated by simple graphs carrying only 5 potentials. 

qH Keywords: quantum graphs, threshold resonance, singular vertex coupling, quantum control 

I 1 PACS: 03.65.-W, 03.65.Nk, 73.63.Nm 

> 

in 

^n 1. Introduction 

in 

^^ Quantum mechanics on graphs is recently attracting more and more attention due to its cur- 

CO rent and prospective use in nanosciences, especially nanoelectronics [1, 2, 3]. In particular, since 

^^ quantum graphs serve as effective models of graph-like structures of submicron sizes, they allow 

^^ to study single electron devices based on interconnected nanoscale wires. One of the first such 

. . applications emerged in the spectral filtering. It is known that a line carrying the ^-interaction is 

usable as a high-pass filter, and similarly, a line carrying the (5' -interaction works as a low-pass 

filter. More complex spectral filters have been developed as well, for instance a trident filter [4], 

;J] or a Y-shaped branching filter, functionning as a high-pass/low-pass junction [5]. 

Cu If a quantum device is easily controllable, its usefulness is markedly enhanced. In principle, 

any device based on a quantum graph, including the aforementionned filters, can be controlled 

by adjusting its internal properties, for example, by varying the parameters of the couplings in its 

vertices. However, such a controllability is hard to be implemented in practice since it requires a 
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real-time modification of a nanoscale object. It is, therefore, highly desirable to develop a device, 
controllable in a simpler manner. 

In this paper we propose a quantum filtering device that can be controlled without afi'ecting 
its internal structure. We consider a model in which constant potentials are applied to the lines 
of a star graph, and the control is achieved through the variation of the potential strengths. The 
design of the device essentially relies on a special type of point interaction used in the star graph 
vertex, namely, on a scale invariant (Fulop-Tsutsui) coupling. It should be emphasized that the 
Fulop-Tsutsui point interaction is very different from usual free coupling and the 6 potential. Its 
exotic scattering properties together with the external potentials on the lines induce a threshold 
resonance effect, on which the filtering function is based. 

It is known that quantum graph vertices with Fulop-Tsutsui couplings can be approximately 
constructed from vertices carrying (^-couplings [6, 7]. In this paper we improve the existing 
approximation schemes, and, in addition, demonstrate its convergence in terms of the scattering 
matrices. Since the vertices carrying the (5-couplings can be approximated by regular, smooth 
potentials [8], the realization of the required Fiilop-Tsutsui couplings is within reach of the 
technology in the near future. Therefore, our quantum graph model of a quantum device is of a 
practical relevance, not just a mathematical abstraction. 

The paper is organized as follows: In Section 2 we prepare the necessary theoretical basis, 
including the scattering formalism of quantum graphs with external potentials on the lines. In 
Section 3 we introduce the concept of a potential-controllable quantum filter. The first device of 
this type is designed in Section 4, namely a band-pass spectral filter with one passband tunable 
by adjusting the potential on the controlling line. Before generalizing the result to more complex 
devices, we devote Section 5 to a practical question. Although the control potential is assumed 
to be constant on the whole line, we demonstrate that its small perturbation in a segment close to 
the vertex, which can be expected in practice, does not distort the filtering function. The idea of 
Section 4 is extended to filters with two controlling lines in Sections 6 and 7. The filter proposed 
in Section 6 has two passbands, and the position of each of them is independently tunable by a 
control potential. The filter developed in Section 7 has one passband of adjustable bandwidth: 
its upper and lower cutoff energies are both tunable by the potentials on the controlling lines. In 
Section 8 we introduce a spectral filter with r controlling lines, r >2. The filter has r passbands, 
each of which can be independently adjusted (as well as disabled) by the potentials on the con- 
trolling lines. This result is achieved by generalizing the idea of Section 6. Finally, a filter with 
multiple independent outputs is studied in Section 9. The outputs work as individually tunable 
band-pass filters; the position of the passband at each output is adjustable by tuning a potential 
on a dedicated controlling line. 

The construction of the Fiilop-Tsutsui couplings is discussed in Sections 10 and 11. In 
Section 10 we demonstrate how they can be approximately represented by a small web carrying 
5-couplings in the vertices and vector potentials on the lines. In Section 11 we calculate the 
on-shell scattering matrix of the approximating graph and prove its convergence to the scattering 
matrix of the required Fiilop-Tsutsui vertex in the small size limit. The paper is concluded by 
Section 12 in which we discuss possible extensions and related open problems. 

2. Preliminaries 

A wave function of a particle confined to a graph having A^ lines Ei, . . . , E^ of the lengths 
^(ei), . . . , {(6^) has A^ components, *!* - {t//i,il/2, . . . , iAat)^, and the corresponding Hilbert space 



'H is given by ^ .^j L?({(Ej)). Let us assume that there are potentials Vi, . . . ,Vn and vector 
potentials Ai, . . . ,An imposed on the graph lines. The Hamiltonian action is given by 



H 



r«Ai 



lAwJ 
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2m 






(1) 



where m is the mass of the particle and q is its charge. Note that if there is no vector potential, 
the Hamiltonian takes a simpler form 



H 



(^i 
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h^ 



f-i^-A'i' + Vi-^i 



h- 



^-t^N + yN-<pA 



(2) 



2.1. Vertex couplings 

Let us consider a graph vertex at which n lines are coupled. We assume that the wave function 
components at these lines are denoted by iffi(xi), . . . , t//n(x„) and that the coordinates are choosen 
such that Xj - corresponds to the vertex for all y = !,...,«. We introduce the vectors 



no) 



(<pm 



MO)) 



and 



T(0) 



(^[(0) 



K(0)) 



(3) 



Since the Hamiltonian is a second-order differential operator, the boundary conditions at the 
vertex couple ^(0) and *P'(0). The most general form of the boundary conditions is 



A*P(0) + BT(0) = , 



(4) 



where A and B are complex nxn matrices. 

To ensure the self-adjointness of the Hamiltonian, which is in physical terms equivalent to 
conservation of the probability current at the vertex, the matrices A and B cannot be arbitrary but 
have to satisfy the requirements 



rank(A|B) = n, 

the matrix AB* is self-adjoint. 



(5) 



where (A\B) denotes the n x 2n matrix with A,B forming the first and the second n columns, 
respectively [9]. The relation (4) together with the requirements (5) describe all possible vertex 
boundary conditions giving rise to a self-adjoint Hamiltonian. 

The requirements (5) can be partly included into the boundary conditions (4) by writing the 
matrices A, B in a special form. One of the ways was discovered by Harmer [10] and indepen- 
dently by Kostrykin and Schrader [11] who transformed (4) with (5) into "compact" boundary 
conditions 

(^( - 7)^(0) + i(^( + I)T(0) = , (6) 



where 'Z/ is a unitary nxn matrix. Other way was proposed in [6]. It consists in transforming (4) 
into the block form 

(^0 )^'^o^"(-r A )^(0), (7) 

where re {0, 1 , . . . , «), /"-"^ is the identity matrix « x n, T is a general complex rxti-r matrix and 
5 is a Hermitian matrix of the order r. We call (7) ST -form of boundary conditions. It is simple 
and unique, but on the other hand it requires an appropriately chosen line numbering, see [6] for 
details. 

In this paper we primarily use the S T-form (7) that best meets our needs. However, any form 
of boundary conditions has its pros and cons and is suitable for certain applications. 

2.2. 6-couplings 

The (^-coupling (or "(J potential") in a vertex of a quantum graph is characterized by boundary 
conditions 

n 

^/0) = ^,(0)=:iA(0), jj=l...,n, ^^;.(0) = a«A(0), (8) 

where a e M is a parameter of the coupling. The (J-coupling is the second most natural singular 
interaction (after the free coupling) by reason of its simple interpretation: It can be understood 
as a limit case of properly scaled smooth potentials [8]. 

Remark 2.1. Strictly speaking, the S potential of the strength Vq iV(x) - V()6(x)) is characterized 
by boundary conditions 

^/0)-^,(0)=:iA(0), jj^l...,n, ^^;.(0) = -^yo«A(0), (9) 



where m is the particle mass. However, since the particle in question is almost eclusively electron 

2m - 



and, therefore, m - m^ - const., the term ff Vq is commonly written as one single constant a. 



usually called parameter. 

2.3. Scattering matrix 

If a quantum particle with mechanical energy E living on the star graph comes in the vertex 
from the ^-th line, it is scattered at the vertex into all the lines. The y-th component of the 
final-state wave function is given by 



¥j\x) = 



Sje^e'^'' forj^ 



(10) 



where Sji are scattering amplitudes, kj are angular wavenumbers on the corresponding lines, and 
coeflicients 1 / -y/fey are involved for proper normalization. For any j, the wavenumber kj is equal 
to 

/ 2m(£ - V,) 

where Vj is the potential on the j-th line. The matrix S - {Sj{] is the scattering matrix of the 
graph. For a normalized wave function coming in from the y'-th line, Sij is interpreted as the 

4 



complex amplitude of transmission into the 7-th line (for j + €), whereas Sa represents the 
complex amplitude of reflection. The matrix S depends, in addition to the internal properties of 
the vertex, on the potentials Vi, y2, . . . , y„, and on the particle energy E. In order to express this 
dependence, we will often use the full symbol 

>S(£;yi,...,y„). (12) 

In case that Y j = for all y = 1, . . . , n, we will simplify the notation by omitting Y j, i.e., 

>S(£)=.S(£;0,...,0). (13) 

In order to derive a formula for S, we define matrices Y and Y' such that Y jt - *? . (0), 
Y'.^ = (^Pf )'(0). With regard to (10), it holds 

F = D ' + D^S = D '(/ + S), 
r = \D^{-D-^ + D-^S) = \D{-1 + S\ 

where D is the diagonal matrix 

D = diag(Vfe^,..., V^), (15) 

for kj given by (11). Any wave function ^(x) = {\p\{x), . . . ,i/'„(x))-^ (the superscript T denotes 
the transposition) on the graph obeys the boundary conditions (4) determining the vertex. In 
particular, b. c. must be satisfied by the final-state wave function T . determined in (10) for all 
j, hence AY + BY' = 0. When we substitute for Y and Y' from (14), we obtain 

AD-\l + S) + iBD(-I + S) ^ , (16) 

which leads to the sought expression for the scattering matrix S(E; Yi, . . . , V„): 

S(E; Vi, . . . , y„) = -(AD-^ + iBD)-\AD-^ - iBD) . (17) 

Matrix S is always unitary, and squared moduli of its elements have the following interpretation: 
\Sj([^ for j i^ { represents the probability of transmission from the /'-th to the j-th line, ISal^ is 
the probability of reflection on the /'-th line. 

Note that if there are no potentials on the graph lines, i.e., Yj - for all y = !,...,«, then 
kj - ^JlmE/h? for all j, and the scattering matrix can be written in the simpler form 



S(E) = -(A + ikB)-\A - ikB) , where k = ^J2mE/^f . (18) 

Equation (18) allows to find a relation between the boundary conditions (4) and their form (6): 
Observation 2.2. The boundary conditions (4) are equivalent to 

(S(Ei)-D^'(0) + i(S(Ei) + DY(0)^0 forEi^^. (19) 

Proof. It suffices to substitute for S{Ei) from (18), and then to multiply the equation by the 
matrix A + iB from left. D 
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2.4. Fillop-Tsutsui couplings 

Vertex couplings with S{E) independent of E are called scale invariant, or Fiilop-Tsutsui, 
couplings [12, 13, 14]. In this paper we use the latter term, since the name "scale invariant" 
sometimes leads to misinterpretations. It can be easily shown that the scattering matrix S(E) 
corresponding to the boundary conditions (7) is energy-independent if and only if the Hermitian 
matrix S in (7) satisfies 5=0, i.e., iff the boundary conditions take the form 



/("^ T \ , / 





j(n-r) 



no). 



(20) 



The scattering matrix corresponding to b.c. (20) can be explicitely expressed using formula (17) 

\ „ //('' T\ 

or 



to which we substitute A 



-T* r(n-r) 
-1 



B 



S{E) 



Un-r) ^ y*tY^ 2t* - (i^"-'^ + r*r) ' (/<""''' - t*t) ^ 



(21) 



This result can be found for example in [7]. 



2.5. Filtering property of the 5-interaction 

The (5-interaction can serve for a simple spectral filtering (Fig. 1). This effect can be demon- 
strated with the help of the scattering matrix. The boundary conditions (8) are equivalent to 



1: input 



6 

-o- 



2: output 



Figure 1 : (5-interaction on tlie line as a liigh-pass spectral filter. 



- (^1 i) **'<^0) + (o o) '^'^'^^ " ^' "^^^^^^ "^^"^ equation (18), 



S(E)^ 



( -ia 2a/£ ^ 

2 Vfi+iff 2 Vl+io- 

2\[E -ia 
V 2V£+ 



2V£+ 



(22) 



where a e M\{0) is the parameter of the (5-interaction. The transmission amplitude input — > out- 
put (1 — » 2) is given by S2i{E), hence we obtain the transmission probability 



P(E) = ISniEf 



4E 



4E + a^ 



(23) 



which satisfies 



P(E) * for £■ «c 0-2 



and 



riE) * 1 for £■ » (3-2 



(24) 



To sum up, the system exhibits fuU reflection for small energies and full transmission for high 
energies, and, therefore, can be used as a high-pass spectral filter. 



3. Potential-controlled quantum device 

The main goal of this paper is to design a quantum fiher that is controllable by an external 
potential, as schematically illustrated in Figure 2. Any such device has, therefore, three (or more) 
lines that are used in the following way: 

• Line 1 is input. Particles of various energies are coming in the device along this line. 

• Line 2 is output. Particles passed through the device are gathered on this line. 

• Line 3 is controlling line. We assume that this line is subjected to a constant (but ad- 
justable) external potential U. Adjustment of the potential directly controls the flow from 
the input to the output. 

The concept can be naturally extended to devices having more than one controlling line, if more 
device parameters are to be controlled. Similarly, devices with several output lines, or several 
input lines, can be considered. We will deal with such cases in Sections 6-9. 



3: control 
U 



1: input 



2: output 



Figure 2: Schematic depiction of a quantum device controllable by an external potential. The transmission through the 
channel 1 — > 2 is being controlled by adjusting a potential U on the controlling line 3. 

It is desirable that the device is technically simple. We satisfy this requirement by designing 
the device modeled by a quantum star graph. In its vertex, we assume a special, Fiilop-Tsutsui 
coupling. The parameters of the coupling will be specified later. They determine the transmission 
characteristics of the device. 

In order to describe the effect of the adjustable control potentials on the scattering properties 
of the graph, we need to generalize the scattering matrix formula (21). 

Proposition 3.1. Let the vertex coupling in the center of a star graph is given by (20) and there 
be constant potentials V\, . . . ,V„ on the graph lines. Then the scattering matrix is given by 

S(E; Vi, . . . , y„) = -/("^ + 2 i J^^<'^. j (e^i) + TQl;r*\' ( 2(1) re(2) ) , (25) 

where 



G(i)=diag|{/l--^,...,Tjl--^ 



e(2)=diag|{/l--^,...,-^l--^ 



(26) 
Proof. The result can be obtained from equation (17) in a similar manner as formula (21). D 



4. Potential-controlled spectral band-pass filter 

The simplest device we are going to construct is based on a star graph with just 3 lines. 
The Fiilop-Tsutsui coupling in its center is described by the boundary condition (20) with two 
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additional assumptions: (i) r = 1, (ii) T is real. The reason for assuming T real will become 
obvious later in Section 10. Hence T - (a b), a,b eR, and the boundary conditions read 



(27) 



The graph is schematically illustrated in Fig. 3. The roles of individual hnes are as stated in 
Section 3:1= input, 2 - output, 3 = controlling line. The scattering properties of the graph 



(1 a b\ 


(^M 




ro 0^ 


(MO)] 





«a;(0) 


= 


-a 1 


«A2(0) 


[o oj 


U3(0)J 




[-b ij 


U3(0)j 



-itofi 



U 



S^^e 



\kX2 



1 2 

FT: (a b) 

Figure 3: Scheme of a quantum spectral filter controllable by an external potential U on the line 3. 

in question are determined by the scattering matrix S{E; 0, 0, U), which can be calculated using 
Proposition 3.1. If we substitute T - (a b) and V\ - V2 - Q,V^ - U into equation (25), we 
easily obtain 



>S(£;0,0,f/) = 



1 



I 1 



\ + a?- + b^^ 



a^-b^^ 


la 


2^V? 


la 


a^-\- b^^ 


2afoV? 


2b yf^ 


labyf^ 


bH-l- 



(28) 



where ^ = yfTHIjE (thus V| = ^1 - U/E). 

A quantum particle with energy E coming in the vertex from the input line 1 is scattered at the 
vertex into all the lines. The corresponding scattering amplitudes are given by the entries in the 
first column of (28). Let us denote for simplicity S .^ (E) = [>S(£';0, 0, t/)]/i for all j - 1,2,3. 
The transmission amplitude input — > output equals 



j(f/) 



S^'^.'iE) 



2fl 



1 +fl2 + 



b^^l 



(29) 



u 

E 



the reflection amplitude is 



:■((/) 



S\\>{E) = 



1 - fl2 _ /,2 f[Z 



l+a^+b^ Jl-f 
and the transmission amplitude 1 — > 3 is given by 

2bf^ 



(30) 



Sf(E) 



l+a^+b^ ^1-f 



©(£ - U) . 



(31) 



The Heaviside step function &(E - U) is added in equation (31) to make the expression valid for 
all energies E, including E < U. It represents asymptotically no transmission to the line 3 below 
the threshold energy Zith = U. 
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We are interested above all in the transmission probability from the input line 1 into the 
output line 2, which we denote by P^^HE). Since P^^\E) = \Sf^\E)\^, we have from (29) 



'P^"\E) = 



4a' 



for E> U, 



(l+a^f+h\U/E-l) ^°^ E < U. 



We observe that for a given constant potential on the line 3, 'P^^\E) as a function of E grows in 
the interval (0, U), decreases in the interval (U, oo), and satisfies 

lim P'-"\E) = , (33) 

lim^(^)(£) = (-^f , (34) 

limP'^^\E)^i ^ ^^ J . (35) 

£^oo yi + a-^ + jyZ I 

Therefore, if the parameters a, b are chosen so that b^ » Aa^, the function P^^\E) has a sharp 
peak at £ = L'^. Furthermore, equation (34) implies that the peak is highest possible (attaining 1) 
for a = 1. To sum up, the device from Figure (3) built for a, b satisfying 

a = 1 and b^ lA » 1 (36) 

has the following property: The transmission probability input — > output is high for particles 
having energies E x U and perfect for E - U, while it is very small for particles with other 
energies. It means that such a device works as an adjustable band-pass spectral filter, allowing 
to control the passband position by the potential put on the controlling line 3. The situation is 
numerically illustrated in Figure 4. 

Let (36) be satisfied. Filtering properties of the device, in particular the sharpness of the peak, 
are determined by the quantity y6 := b^ 12. With regard to the above results, it holds P^^XE) » 
for £ ^ and P'^'^\E) » 1/(1 + f^f for £ ^ oo. The bandwidth W of the filter, i.e., the width 
of the interval of energies E for which !P*^'(£') > 1 /2, can be calculated from equation (32). The 
bandwidth depends on U and/?^, and equals 

2(2- V2)yS2 

W=- !^ ^ U. (37) 

(y32-3+2V2)062 + l) 

Since yS^ » 1, it holds W ~ 2(2 - V2)f//j6^ ~ l.nU//3^. 

Remark 4.1. The resonance at the threshold energy Eih = U is related to the pole of the scattering 
matrix which is located on the positive real axis at 

b* 

M - (1 +a^y 

on the unphysical Riemann surface which is connected to the physical Riemann surface atE = U. 
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Figure 4: Scattering characteristics of the graph from Fig. 3 with parameters a = \, b = A (i.e., /? = 8), plotted for 
the control potential set to J/ = 1. The transmission probabihty P^^^E) as a function of E is plotted in the left figure. 
The right figure shows the reflection probability \S\j (E)\~ and the probability of transmission to the controlling line 

\S\"\Ef. 



Remark 4.2. A band-stop spectral filter can be constructed in a similar way. If we begin with 
r - 2 and 7^ = , I and choose c,d eR such that c - d:^ 1, we obtain 



P^"\E) = 



4c4 



2 ' 



(39) 



hence 'P''''(£') ~ 1 every where except a certain narrow interval around £ = f/, andlim£^j/!P*^'(£') 
0. 



5. The role of the control potential 

The filtering property of the device proposed in Section 4 has been derived on the assump- 
tion that the control potential U is constant on the whole line 3. A natural question arises: What 
happens when the control potential is not constant precisely up to the junction and/or U{x) is 
supported only by a certain segment of the line 3? In this section we show that under certain 
conditions, the filtering property is not distorted. To obtain the result, we prove, for a general 
potential U{x), that the transmission probability p(^''^^\E) can be expressed in terms of the re- 
flection amplitude calculated for the potential Uix) on the real line. The formula we obtain will 
allow not only to answer the above question, but also to better understand the function of the 
filter. 

In order to find the transmission amplitude S'^^ \E) for a given potential U{x) on the 
controlling line, we need to apply boundary conditions (27) on the wave function components 
i/'i(x), t//2ix), il/T,{x). The particle motion on the lines 1 and 2 is free, thus the boundary values 
1/^1(0), i/'j(0), (/'2(0), (^2(0) at the vertex can be expressed in terms of the scattering amplitudes, cf. 
equation (10), 

'Ai(0) = -^(-1+^11), ^2(0) 



i|f\{0)^i^fk{-\+Sn), 



Tk^^' 



^'(0) = iV^^21, 



(40) 



where S\\ is the reflection amplitude to the input line 1, >S2i is the transmission amplitude to the 
output line 2, and 

k-^. (41) 
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Our next aim is to express the boundary values i/'3(0) and ifr'^(Q) in a convenient way. For that 
purpose, let us consider a line with the potential 



U(x)^ 



U{x) for jc > , 
f or a; < , 



(42) 



see Figure 5. For i/'3(x) being the given wave function component on the controlling line, we 

define the function 

{iffiix) for X > 0, 

cUe'*- + iR^e-'*-) forx<0, ^^^^ 

where C and K are determined by the continuity of 4r(x) and tfr'{x) at jc = 0: 



^c(i+??) = iA3(0), iV^c(i-??)-^^(0). 



(44) 



The function ijr(x) is obviously a solution of the Schrodinger equation for a particle on the 



C(e + Re VVk '/'aW 



U{x) 



Figure 5 : A line with tlie potential U(x) in the part .v > 0. 

line with the potential Uix). The term C(e''^ '/ V^ + "Re"'*'/ V^) corresponding to x < in 
definition (43) implies that ij/ix) is the wave function of a particle of energy E moving along the 
line from the part x < to the part x > 0, and !R represents the reflection amplitude at the point 
x = 0. 

Remark 5.1. From equations (44) it foflows "R = [ik>J/3(0) - tfr'^(0)) / (iki//i(0) + iJ/'^iO)). This 
formula can be generalized. It is easy to show that for any xq > 0, it holds 



nxo) 



ikil^3(0) - 1/^3(0) 



(45) 



where 'R(xo) represents the reflection amplitude corresponding to the shifted and truncated po- 
tential Uxoix) - U(x - xq) ■ @{x - xq), measured at the point x = 0. 

Now we use equations (40) and (44) to subsitute for (/'>(0) and (A'(0), j = 1,2,3, in the 
boundary conditions (27). After an easy manipulation we obtain the system of equations 



(1 a b(\-K) 
fl -1 

b -{l+K)} 



•S21 

c ) 



( l" 

a 
b] 



(46) 



therefore, the transmission amplitude input — > output equals 



'21 



(£) 



20(1+-^) 



{\+a^)(\+'R) + bH\-'R) 
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(47) 



If a = 1 and/3 :-b^ 12 as introduced in Section 4, then 



l+'R+Pd-'R) ' ^"^^^ 

Equation (48) relates the effect of the control potential, encapsulated in the reflection amplitude 
'R, and the actual transmission probability y5(^(^))(£') of the filter. 

Since /? » 1 by assumption, we observe from equation (48) that the filter is open for the 
particles that would be reflected by the potential U(x) with the amplitude K ~ +1, and closed for 
particles with ^(1 -K)] » 1. We can say that the control potential opens the filter by reflecting 
the particles entering the line 3 back to the vertex with the amplitude +1. In particular, if the 
particle has high enough energy to pass through the potential barrier U{x), we have K ~ 0, hence 

p(UM)(E) a: 1/(1 +j6)2~0. 

In case that the potential is finitely supported and/or separated from the vertex by a gap, the 
reflection amplitude K is usually oscillating, thus ^(^(^W is oscillating as well. Furthermore, if 
the constant part of the potential is long enough, the oscillations are very rapid and, therefore, 
imperceptible (the "semiclassical" short wavelength limit). In such cases we mollify the quickly 
oscillating function ^(f'('))(£) by a convolution with the Friedrichs' moUifier w^ for a certain 
small e > 0, where 

w,(z) = | Kie"*" fork|<e, (49^ 

I otherwise 

r\ 1 

and Ci is chosen such that J Cie '-^^ dz = 1. As we will see in examples below, the obtained 
mollified function 

^^^«>(£) = ^*^(^W))(£) (50) 

is essentially identical with the transmission probability input — > output in the ideal case when 
the control potential is constant on the whole line 3. 

Example 1: Control potential supported by a finite segment [0, L\ 

If t/[o,L] - U ■ X[0,L] ioi xi(),L] being the characteristic function of the interval [0,L], we have 

^ (l-^sin^fcL 

^"•^^ ( 1 + ^2) s in ^1^1 + 2i^ cos ^/tL ' 

where ^ = VI ~ U/E. It can be shown that the transmission probability, given by equation (48) 
with K obeying (51), grows in the interval of energies [0, L'^] . Furthermore, if the support of L'^(x) 
is long enough in the sense L :s> hi yjlmU, then kL :^ 1 for all particles with energies E > U, 
and one can prove that 

• lim£^[, 'P(^M>(£') » 1, 

• in the region of energies E e [U, 00), the transmission probability f'(^[»ii)(£'^ quickly oscil- 
lates between P^i^ » 1/(1 +/3)2 » and !Pn,ax ~ 1/(1 + (1 - U/E)0. 

The situation is illustrated in Figure 6 (left). The figure shows also the mollified transmission 
probability iP^^Jj^'^E) = (w, * 'P''""'»^)(E). The function 'P^^T^^^') is very close to the transmis- 
sion characteristics obtained in Section 4 for the control potential constant on the whole line 3. 

The condition L » fi/ ^j2mU = An/ (In) means that the support of the control potential 
is substantially longer than the de Broglie wavelength of the particle with energy U. Such a 
condition is usually easily satisfied by any macroscopical length L. Indeed, if we consider an 
electron and, for instance, U - 0.1 eV, then h/ y/lmU x 2 • lO"'" m = 0.6 nm. 
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0.0 Jf 




Figure 6: The transmission characteristics of the quantum filter designed in Sect. 4, controlled by a potential supported 
by the finite segment [0, L] (left) and [f, L] (right). The mollified transmission probabilities are plotted by thick black 
lines. The parameters are chosen as J/ = 1, L = 2000, { = 1. If we assume the particle to be an electron, the energy 
scale to be given by meV, milielectronvolt, and the length scale to be specified by nm, nanometer, then J/ = 1 meV 
corresponds to the wavelength A = h/ V2»J(. U = 38.8 nm. For L 2> A/(27r) 2> d, the mollified transmission probability is 
essentially identical to the transmission characteristics obtained in the ideal case L = oo, cf. Fig. 4. 



Example 2: Control potential supported by a finite segment [£, L], ( > 

In the second example we focus on a control potential with the support separated from the 
vertex by a gap of the length d, i.e., t/[f,L] - U ■ X[C,l]', the corresponding coefficient "Ricl] equals 



%e.u - 



(1 -^^)sm^k{L-0 



'■^ ( 1 + ^2) sin ^^(^ _f) + 2i^ cos ^k(L - () 
If the lengths L, € satisfy 



^2ikt 



L» 



V2m[/ 



and 



h 



<K 



V2^' 



(52) 



(53) 



then the values ofR^ci^ are close to the values ofR^Qi^^ obtained in Example 1 (cf. eq. (51)). 
Indeed, 



for E < U oi E ^ U,we have k 
equation (53), hence e^'*^^ ^ 1 ; 



y/lmE/h < y/lmU/h, hence 2k£ <s 1 due to the second 



for £■ » f/, it holds ^ a; 1, hence "R^f ij * ~ "R, 



HO.i]- 



Therefore, the transmission characteristics of the filter with U[[^l] - U ■ X{t,L] are similar to those 
described in Example 1, see Figure 6 (right). Namely, the transmission probability 'P^^iw\E) 
is quickly oscillating, but the mollified transmission probability P^^'^^(E) - (co^ * !P*^i'^'-i')(£') 
well approximates the ideal characteristics obtained in Section 4. The conditions (53) mean that 
{ <c Au/i2n) <K L, where Au - h/ y/lmU is the de Broglie wavelength of the particle with 
energy U. 

In practical implementations, the profile of the control potential U{x) might differ from those 
studied in Examples 1 and 2 above. In particular, the potential may be constant in the interval 
[£, L] and taper off to zero in the segment [0, £]. Generally, if the conditions (53) are satisfied, 
the result is essentially independent of the shape of the fall-off of U{x). To demonstrate this, let 
us distinguish two situations. 
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If £ < U or E ^ U, the local wavelength A(x) - hj ■\J2m{E - U{x)) of the particle in every 
point X e [0, {] satisfies A(x)/(2n) » { (the "long wavelength limit"), hence i/'3(0) ~ i/'3(^), 
i/'jCO) * i/fjC^). Therefore, with regard to Remark 5.1, 



n^ 






= 5R| 



[Q,L-e] , 



(54) 



where ^[o,l-^] is the reflection amplitude corresponding to the potential t/[o,L-^] = U ■ 
X[Q,L-e] ■ Since moreover L-{ » Au/(2n), it follows from Example 1 that the corresponding 
mollified transmission probability is similar to the transmission characteristics found in 
Section 4. 

• If £■ » U, the particle is almost fully transmitted through the potential barrier, hence 
"R a 0. Therefore, y5(f^W)(£') » 1/(1 +^)2 due to equation (48), which coincides with the 
result obtained in Section 4. 

To sum up, the filter works well also for a finitely supported potential which is not constant 
precisely up to the junction, on condition that the support of the potential is long and the segment 
in which the potential tapers off is short, both with respect to the de Broglie wavelength of the 
particle with energy U. A similar result holds true also for the other filtering devices described 
in subsequent sections. 



6. Spectral filter with two passbands 

In this section we design a band-pass quantum filter with two passbands such that their posi- 
tions are independently controllable by two external potentials. 

The device will be naturally based on a star graph with 4 lines that have the following mean- 
ing: 1 - input, 2 - output, 3 and 4 - controlling lines subjected to constant external potentials 
U, V; see Figure 7. 

The vertex coupling in the graph center is of the Fiilop-Tsutsui type in accordance with the 



concept introduced in Section 3. This time we start with an ansatz r - 2, i.e., T = 



as in the previous section, assume that T is real, thus a, b,c,d e 
described by the boundary conditions 



a 



I, and. 



1 





a b\ 


(^\{% 





1 


c d 


«a;(0) 











iA^(O) 








oj 


W;(o)i 



CO 











-a 


-c 


\-b 


-d 






0^ 


i^m\ 








«A2(0) 


1 





<A3(0) 





\) 


U4(0)J 



c dj 
Therefore, the vertex is 



(55) 



W,V). 



The calculation of the transmission amplitude vSjj' (E) = [S{E; 0, 0, U, y)]2i can be per 



formed in the same way as in Section 4 using formula (25); we find 

-2(«c^/r^+M^r^) 



sfr\E) = 



1 -H (fl2 -H c2) ^l-^+(b^+ d^) yjl-l+iad- bc)^ 7^ - f ^ 



(56) 



The transmission probability input — > output is given as !P*^'^'(£') = \S^^' (E)\^. Since the filter 
we seek shall have two passbands at £ = t/ and E - V,we require zero transmission probability 
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Figure 7: Scheme of a quantum spectral filter controllable by two external potentials U and V on the controlling lines 3 
and 4, respectively. 



for £■ — > 00 and high transmission probabiUties (preferably » 1) at £ = U and E = V. The first 



requirement (i.e., hm^^oo ^Sjj' (E) = 0) leads to the condition 

ac + bd -0. 



(57) 



Then we proceed to the requirement 'P^^-^\U) ~ 1, 'P'^'^\V) x 1. Let us assume without loss of 
generality V < U. Since it holds 






2bd 



f^v ] 



4oV(^-l) 

l+(fl2+c2)2(£-l)' 



we arrive at the condition 



fl^ = c^ » 1 . 



^2 = £/2 » 1 , 
Conditions (57) and (60) together lead to eight possible expressions for T, namely 



a a \ la a\ la -a\ I -a a 
a -a]' ~\-a a/' ~\a a ]' ~\ a a 



for a > such that a^ ^ I. We choose the first of them, i.e.. 



T = 



(58) 



(59) 



(60) 



(61) 



(62) 



but it is easy to show that all the solutions (61) result in the same transmission probabilities, thus 
a different choice would not make any difference. Now we use the formula (25) to calculate 
the corresponding input — > output transmission amplitude S^i' (E), the reflection amplitude 
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S.' '(£), and the remaining transmission amplitudes .si ' (E),S.' (E): 



ST(E)- 






ST(E) 



ST(E) 



ST(E) 



1-4«^V^V^ 



(i+2«2^rr|)(i+2fl2^/rr|) 



^ 



2fl*l 



1 + 






: ®(E - U) , 



la^^fT^ 



@(E - V) . 



1 + 



The input — > output transmission probability is equal to P^^'^\E) - \S\' (£)P, thus 



p'^'^Ke) 



(l+4«4(g-i))(i+4«4(|-l)) 

4fl^ U-V 

(l+4fl4(y_i))|i+2a2-'^^''' -^ 

2a2(Vr7-' 



l+2a2 -v/l-¥ l+2a2 V^T 
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for £■ < y, 
for y < £■ < f/, 

for £■>{/. 



If t/, y satisfy 20^ yfT^VjU » 1 and 4fl'*(f//y - 1) » 1, then 

lim p(^-^>(£) = ^ ^ ;^ 1 , 

E-^y l+4fl4(^-l) 

Um r^^'^^E) = ^ ^ * 1 , 

hm 'P^"'^\E) = , lim P^"'^\E) = , 



<p(f/,v) (-£■-) ^ Q £qj. ^jj £■ except for certain small neighborhoods of U and y. 



(63) 



(64) 



(65) 



(66) 



(67) 



(68) 

(69) 

(70) 
(71) 



Therefore, the function P^^-^XE) has two sharp peaks, located at the energies E - U and E = V, 
see Figure 8 (left). 
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The situation is different when one of the potentials vanishes. If V = and U > Q, then 

1 



= 0, 



hm 'P^"-^\E) = 
Hm !P<^'^>(£) = , 

lim r^^'^HE) = ^"' , ^ 1 , 
£->E^ (l+2fl2)^ 

P^^'^\E) ^ for all £ except for a certain small neighborhood of U. 



(72) 
(73) 

(74) 
(75) 



Therefore, turning one control potential down to zero (e.g., V = 0) does not result in a peak at 
£■ = 0, but causes the peak corresponding to V to vanish, whereas the peak located at £ = L'^ 
remains essentially unaffected; see Figure 8 for an illustration. 

To sum up, the device constructed according to the scheme in Figure 7 for T given by (62) 
with 2a^ » 1 works as a band-pass spectral filter with 2 passbands. The positions of the pass- 
bands are directly controlled by the external potentials on the controlling lines, and moreover, 
one band can be suppressed by setting one of the potentials to zero. 





I5iil^ 0-5- fth2=l/f 



1531 1 0-5- 
0.0-r 

\SJ 0-5- 

0.0^ 



i= 



r" 



a = 4 



J^ 



fthl = u 



1.0 

E 



I5„l 0.5 
0.0 

153/ 0-5 
0.0 

154/ 0-5 

0.0-1: 



r^ 



EiU = U ••, 



^ 



1.0 
E 



a = 4 



Figure 8: Characteri.stics of the spectral filter obtained from the graph in Fig. 7 for a = fo = f = -rf = 4. The left figure, 
plotted for U = 1, V = 0.5, illustrates the standard dual-band regime. The right figure, plotted for U = 1, V = 0, shows 
the eft'ect of setting one of the potentials to zero; the filter is switched to the single-band regime. The top part of every 
figure displays the transmission probability 'P^^-^\E) as a function of E, the lower part shows the reflection probability 
|S . ,' (E)\^ and the probabilities of transmission to the controlling lines |S i.' (E)\^ and \S ./ (£)p. 



It can be shown that the bandwidth of the passband at the energy [/ > satisfies Wu ~ 
(1 - 1/ V2)t//fl'^. Similarly, if V > 0, then Wi/ = (1 - 1/ V2)y/fl''. 

The peaks of P^^'^^E) at /ithi = U and Efi,2 = V are related to the poles of the scattering 

"''° ~ 4aT^ ^' respectively. 



matrix in the unphysical Riemann plane at /i^" - ^^ U and E^"' - 
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7. Fully tunable band-pass spectral filter 

The filter designed in Section 4 has one controllable parameter, namely the passband position 
which is adjustable by an external potential put on the controlling line. The aim of this section is 
to construct a filtering device that allows to control both the passband position and the bandwidth. 

Similarly to Section 6, the goal will be achieved using the graph depicted in Figure 7 with 
the Fiilop-Tsutsui vertex coupling given by boundary conditions (55). The control potentials U 
and V will adjust the upper and the lower cutoff energies. 

The transmission amplitude input — > output follows from the results of the previous section, 
and is given by (56). Since the sought filter shall be of the band-pass type, we require zero trans- 
mission probability for £ — > oo, i.e., lim^^oo S^^ (E) - 0. Hence we obtain the condition (57), 
i.e., ac + bd - 0. 

Now let us consider special situation where the lower cutoff energy is zero; in this regime the 
device shall work as a low-pass spectral filter. If we set V = in (56) and use (57) to simplify 
the numerator, we obtain 



\+b^ + (f + {a^ + c^ + {ad - bc)^) J]~^ 



4f '(£) = ^ ■ '- ^= . (76) 



Note that for E e [V, U] = [0, U], the expression VI - U/E is imaginary, and consequently, the 
transmission probability for £ < L'^ (in the intended passband) equals 

^(^•«>(£) = \S^^;'\E)\^ = 5^ 5^ ^ . (77) 

il+b^ + c/2)2 + (fl2 + c2 -H (ad - bcf)^ (l - f ) 

We observe that there is a special choice of a, b, c, d that leads to a constant transmission prob- 
ability in the whole interval (0, U), i.e., to a flat passband. Indeed, if the parameters a, b, c, d 
satisfy condition (57) and at the same time 

l+b^+d^ ^a^ + c^ + (ad- bcf , (78) 

then 

,„o, Ab^d^ I 2bd \^ 

'P^'^E) = = for all £ < f/ . (79) 

(l+b^ + d^y \l+b^ + d^l 

Therefore, we impose both conditions (57) and (78), and require that the transmission probabil- 
ity (79) is as high as possible (our aim is to minimize the attenuation for E inside the passband). 
In other words, we are to solve the optimization problem 

/ 2bd f 

maximize F(a, b, c,d) -\ z r subject to eq. (57) & (78) . (80) 

\ 1 +b^ + d^ j 

At first, let us simplify the problem by eliminating a and c. We express c from (57), substitute it 
to (78), 

\+b^ +d^ ^a^ + {bdlaf + {ad + b^djaf , (81) 

then multiply this equation by a^ and rewrite it as a polynomial equation in a: 

(1 + d^)a'* + (2b^d^ -\-b^- d^)a^ + b^d\\ +b^)^0. (82) 
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This is a biquadratic equation for a with parameters b, d; it can be shown that it has a real solution 
a if and only if b, d satisfy 

l+b^ + d^ - 8b^d^ > . (83) 



Therefore, the optimization problem (80) is equivalent to the problem 

\2 



maximize F{b, d) = 



2bd 



\+b^ + d^ 
One can easily find its solution: 



subject to \ +b^ + d - ^b^d > 



= |d/| = l/V2. 



(84) 



(85) 



We substitute b, d from (85) into conditions (57) and (78), and in this way find a, c. Altogether, 
we obtain eight solutions for T, namely 



±ll 1 
"V2I1 -1 
We choose the first of them. 



V2(-l l)'"V2(l 1 



±^l , ,1. ±^L " I. ± 



1 1 
V2\l 1 



(86) 



(87) 



V2I1 -1; 

but the choice actually makes no difference, as all the solutions (86) lead to the same transmission 
probabilities. 

We observe that the solution (87) is nothing but a special case of the matrix T found in 
Section 6 (see eq. (62)) for a = 1 / V2. Consequently, the corresponding input -^ output trans- 
mission amplitude iSjj' (E), the reflection amplitude S\^ (E), and the remaining transmission 
amplitudes Si.' (E), S\,' (E) immediately follow from equations (63)-(66), i.e.. 



ST(E) = - 



(l+V^)(l+V^)' 



(88) 



etc. Let us assume without loss of generality that V < U. The input — > output transmission 
probability, given by P^"-^\E) = \Sf{^\E)\^, follows from equation (67) for a = 1/ V2: 



j>iu,V)^E) 






/Wi-Ijii+Vi-I;y 



EL 

uv 



for E<V, 
foiV <E <U, 

for E>U. 



(89) 



The function P^^'^XE) has the following properties: 
• If y = 0, then 

P^ufi^E) = 1/4 for £ e (0, U) , 

P^^'^\E) as a function of E quickly falls off to zero at £ > U , 



lim !P(^'°'(£) = . 



(90) 
(91) 
(92) 
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IfO< y < C/, then 



y 



'P^v.v)^V) = 1 - - , 'P^''-^\U) = 



(-^) 



r^"'^\E) decreases in (y f/) , 

'P^^'^\E) quickly decreases for E > U and grows for E < V , 
lim r''"-^\E) = , lim P'"'^\E) = . 



£■— >co 



The behaviour of ^'^^'^^(E) is illustrated in Figure 9 in both situations y > 0, V = 0. 



(93) 

(94) 
(95) 
(96) 
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Figure 9: Characteristics of the band-pass spectral filter obtained from the graph in Fig. 7 for a = h = c = -d = 1/ y2. 
The left figure illustrates the situation U = \ and V = 0.5, the right figure the situation U = \ and V = (the flat- 
passband low-pass regime). The top part of every figure displays the transmission probability ?'"'•'''(£) as a function 
of E, the lower part shows the reflection probability \S\ ^ (£)P and the probabilities of transmission to the controlling 



;(f.V)^E-M2, 



lines|S!,V"(£rand|S^';'''(£)| 



Remark 1 A. The term "passband" is used in this section in a weakened sense, because the trans- 
mission probabiHty for E e [y U~\ mostly does not exceed 1/2. However, considering the char- 
acteristics of the filter (cf. Fig. 9), it is appai'ently reasonable to regard the interval =s [y [/] as a 
passband, especially in the regime V -Q. 

Remark 7.2. If the extremality requirement (79) is left out, conditions (57) & (78), inducing the 
flat passband in the low-pass regime V -Q, have various other solutions, such as 



T = 



a 


Vl-fl2) 


yl\-a^ 


-a ^ 


1 / a 


a \ 


^\\la 


-lla] ' 



for fl 6(0,1) => !P(^'"*(£)| 
for fl > => 'P<^'"'(£)| 



£e(0,(/) 



^a\l-a^). 



Ee(0,U) (^2 + 1)4 



(97) 
(98) 
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Potential-controlled quantum sluice-gate 

As we have seen, when there is no potential on the line A (V - 0), the device behaves 
as a low-pass spectral filter with a flat passband that transmits (with the probability of 1/4) 
quantum particles with energies E e [0, L'^] to the output, whereas particles with higher energies 
are diverted to the other lines, mainly to 3 and 4. This property allows to use the device also as 
a quantum flux controller: If many particles described by the energy distribution p{E) are sent 
along the line 1, the flux J to the line 2 is given by 



X' 



7(t/)= p{E)V'''"\E)AE . (99) 

Assuming the Fermi distribution with Fermi energy Ep larger than our range of operation of 
U, we can set p{E) - p - const. With the approximation 'P^'^''^\E) » i0(f/ - £), we obtain 
J{U) » ^pU, which indicates a linear flux control. Therefore, the device in the operation mode 
y = can be used as a quantum sluice-gate, linearly adjustable by the potential U applied to the 
controlling line 3. In this regime the line 4 serves as a drain. 

8. Spectral filter with multiple passbands 

In this section we generalize the idea of Section 6. For any r > 1 we design a spectral 
filter which has r passbands such that their positions are directly controllable by potentials on r 
controlling lines. 

The design is based on a star graph with n -2r lines, as depicted in Figure 10. The individual 
lines have the following meanings: 

• Line 1 is input. 

• Line 2 is output. 

• Lines 3, . . . , r are auxiliary lines; their possible use will be discussed later on. 

• Lines r+ 1 ,...,« are controlling lines, subjected to adjustable external potentials U\, . . . ,Ur. 

The Fulop-Tsutsui vertex coupling in the graph center is given by the boundary conditions 
written in the S T-form (20), and the matrix T is chosen as T = aG for an a > and a certain 
unitary matrix G or the order r, i.e.. 



For any j,{ < r, j 4^ I, the transmission amplitude from the /'-th line to the y'-th line 

,(t/) 
it 



can be easily obtained from equation (25). If we denote, for the sake of brevity, S ■. (E) = 



[S(E;0, ... ,0,Uu ■■■ ,Ur)]je, then 



je 



^(;^)(£) . ^SnSn ^ 2gj2ga ^_^ 2gj,g,, ^^^^^ 

1 + fl2 ^1-^ 1 + fl2 ^1- f 1 -t- fl2 ^1- f 
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Figure 10: Scheme of a quantum spectral filter based on a n = 2r star graph. The transmission from the intput 1 to the 
output 2 is controlled by r external potentials f/| ,...,{/,- on the controlling lines r + 1 , . . . , n, respectively. 



The transmission amplitude input — > output corresponds to the choice { - 1, j = 2 in (101). It 
holds 



AU)t 



lim S\{{E) = 2^21^11 + 2^22^12 + ■ 



■2^2r^ = 2[GG*]21=2[/W]2l=0, 

and if the control potentials f/i , . . . , f/^ are all nonzero and mutually different, then 



Um S^"\e) = + + ---+0 = 0, 



Iim4^'(£) = 2g2;^+;^ 



'^giegu 



\<(<r 1 + fl2 



aFI 



(102) 

(103) 
(104) 



If fl » 1 and, moreover, the "isolated potentials" condition 

»1 forci\\j,e^\,...,r,j + €. 



Ui 



is satisfied, we obtain 



and at the same time 



limSf^\E)^2g2jgT] forall7 = l, 



E^Ui 



S21 (E) X for all E except for certain small neighborhoods of Ui ... ,Ur 



(105) 



(106) 



(107) 



To sum up, the absolute value of the function iSjj (£) for a » 1 has sharp peaks at £ = 
Ui,...,U„. Since the input -^ output transmission probability equals I^Sjj (ii)!^, the studied 
device works as a spectral filter with multiple passbands positioned at the energies determined 
by the control potentials. The situation is illustrated in Figure 11 in the case of the device con- 
structed for 

n 1 1 i\ 



1 



1-11-1 
11-1-1 
1-1-1 I ) 



(108) 



The peaks of !P*'^'*(£') at the energies E - Uj (j - I, . . . , r) are related to the poles of the 
scattering matrix in the unphysical Riemann plane at E^.° - -^^—^ Uj. 
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Figure 11: Characteristics of the filter with multiple passbands, obtained from the graph in Fig. 10 with the parameters 
given by eq. (108). The control potentials are set to U[ = 0.1, f/2 = 0.25, f/3 = 0.5, U4 = 1. The left figure shows 
the transmission probabilities \S . (E)\^ between the lines 1 — > ;'; due to the choice of G, all the peak heights are 

almost identical. The right figure shows the reflection probability \S\^ (E)\^ and the probabilities of transmission to the 
controlling lines \S[. (£)P for £ = r + 1,. ..,n. 



On the heights of the peaks and the choice ofG 

Let Ui, . . . ,Urhe all nonzero and satisfying the "isolated potentials" condition (105). With 
regard to equation (106), the heights of the probability peaks at £ = Uj equal 



P^^^'^HUj) = |4^'(f/,0l' * 41^2/1^1/ : 



(109) 



thus, in general, they may be mutually different. However, if all the entries of G have the same 
moduli (equal to 1 / ^fr), then all the peaks are (approximately) of the uniform height 



P''^'>iUj) ^ -. . 



(110) 



A unitary matrix G with this property is nothing but G - 4=//*'^', where H^''^ is an Hadamard 
matrix of the order r. Regarding the existence of Hadamard matrices of order r, it holds: 

• a complex Hadamard matrix exists for every r e N, for example the matrix with the entries 

[//(%, = e"0-W-i), (jj^) 

• a real Hadamard matrix (having elements ±1) of the order r is conjectured to exist if and 
only if r = 1, r = 2 or r is a multiple of 4. 

Consequently, a filter with r passbands of the equal height 4/r^ can be constructed for every 
r > 2, but for most numbers r the corresponding matrices G need to be complex. 

Since the studied device is a generalization of the filter designed in Section 6, it has similar 
properties. A more detailed analysis follows. 
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Reduction of the number ofpassbands. The studied filter allows not only to control the passband 
positions, but also to reduce their number to any r' < r, which is achieved by setting r-r' control 
potentials to zero. Indeed, if r-r' controlling lines carry zero potentials (we may assume without 
loss of generality Ur'+\ = t/r'+2 = ■ ■ ■ = t/r = 0), then 



Sf,\E) . X 



^82{8i 



1 



t=l 1 + fl2 ^ 



E 



l+fl2 



Z ^Siegir « X 



^82(81 



e=r' + l 



e=i 1 + fl2 ^ 



(112) 



E 



with regard to the assumption a » 1. Hence, the effect of setting Ur'+i - Ur'+i - ■ ■ ■ - Ur - 
essentially consists in suppressing r-r' peaks. More precisely speaking, !P'^'^"(0) is very slightly 
increased from to 

4 



E^O 



(E) 



(1 + fl2)2 



and the heights of the remaining r' peaks at Uj, j 

^%2ll'l^l'tO 



{=1 
1,.. 



(113) 



lim !P(I^'>(£) ~ 4 

E^Ui 



82j8ij + 



l+a^Z 



r', are very slightly changed from 

2 

82e8ie > (114) 



but the effect is not significant due to a » 1 . Note, however, that reducing the number of pass- 
bands in this way (when the device is in operation) is not equivalent to using a device constructed 
for a smaller r, because the peak heights generally depend on r (cf equation (1 10), related to the 
Hadamard case, and the discussion above it). 

If all the control potentials are zero, then ;P'^''^"(£) = 0, which follows from equation (101) 
together with the unitarity of G. 

Effect of identical potentials. When several control potentials are set to the same value (or if 
their values are very close), it results in the merging of the corresponding passbands into one 
single passband, and, therefore, to the reduction of the peaks of 'P*'^"(£'). In contrast to the 
reduction achieved via setting several potentials to zero (see above), this manner influences the 
height of the merged peak. If, for instance, U\ - ■ ■ ■ - Uj (j < r), then 
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therefore, the height of the peak at £ = L'^i is changed from » 4|g2i^iil (see (109)) to 

2 



f'<'^"(t/i) ~ 4 



z 

(=1 



82t8u 



(116) 



which may be a significantly different value, depending on the matrix G. The other peaks remain 
essentially unaffected. 

If all the control potentials are equal, it holds 'P^^'^^\E) = due to equation (101) and the 
unitarity of G. 

A similar situation happens when some of the potentials are close to each other, U\ ^ . . . ^ 
Uj (j < r). 
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Remark on the auxiliary lines. It follows from equation (101) that the lines 3, . . . , r can in prin- 
ciple serve as output lines. If the particles passed throught the vertex are gathered on all of the 
lines 2, 3, . . . , r, then the total (aggregated) transmission probability is 



therefore, the cummulative peaks are higher, namely 



(117) 



limnT(£)-4;^|g2/|gi,-|- 

J=2 



E^Ur 



(118) 



For instance, the cummulative peaks in the Hadamard case G - (1/ y[r)H^''^ satisfy !Pj^j (Uc) ~ 
4(r - l)/r^. If r » 1, we have !P,q( (Ue) ~ 4/r, which is a significantly higher value than the 
single-channel peak heights 4/r- found in equation (110). 

9. Spectral branching filter 

In this section we study a filtering device with multiple independently controllable outputs. 
The device is based on a star graph with n - 2r lines, see Figure 12: 

• Line 1 is input. 

• Lines 2, . . . , r are outputs. 

• Line r -l- 1 is drain. 

• Lines r+2, . . . , n are controlling lines, subjected to adjustable external potentials U2, . . ■, Ur. 




^ a O-a 



Figure 12: Scheme of a quantum branching spectral filter based on a n = 2r star graph (for r = 4). For every j = 2,. . . ,r, 
the transmission to the output line ;' is controlled by the external potential U j on the line r + j. 

In the center of the graph there is a Fiilop-Tsutsui vertex coupling given by the boundary condi- 
tions 
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where a > is a parameter such that a^ » 1 . 

For any j - 2, . . . , r, the transmission ampHtude of the channel 1 — > j, denoted by S ^ (E) = 
[S(E;0,...,0,Uu...,Ur)]ji,is 



Sf(E) 



(l+flV)(J,+^,.) 



( 1 + 1:^=2^ 



r 1 \ 



+ff 



1 + y -^ 



+ff 



(120) 



where ^i— a/1 - -^ for all ^ = 2, . . . , r, and the transmission amplitude input — > drain equals 

(121) 



z(V) 



2a 



With regard to equation (120), it holds 

lim S'-^\E) = V; = 2, . . . , r. 

Let the nonzero control potentials be mutually different in the sense 

» 1 for all j, ^ = 1, . . . , r, j ?i € such that Uj + 0, 



1- 






(122) 



(123) 



and let /iq be the number of the control potentials set to 0, i.e., /zq = #{^ e {2, . . . , r) | f/^ = 0). 
Then, with regard to the assumption a^ » 1, 



lim5<,^'(£) = 



1 + /zo + 4 



£^0 ^'^ 1 + a^r \ +k+ '^^~^° 
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n,^^ = 



-2 -2 



+ r A- 
and furthermore, for every j, I such that U j + ^ + U(, it holds 
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The transmission probability to the output line j is given by P ■ (E) = \S .^ (E)\ . If the condi 
tion (123) is satisfied, then for every j - 2, . . . ,r such that Uj + 0, we have 



lim !Pf '(£) = , 
lim 'P^"\E) ^ , 

lim 'Pf (£) ;* 4 ' 

lim 'P^^\E) ^ ^ 



for {/^ ^ , 



(127) 
(128) 

(129) 
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and P .. \E) » for all E except for a certain small neighborhoods of nonzero potentials U( + 0. 



ji 



3(f^)/ 



In other words, if U j + 0, then the probability V (E) has a principal sharp peak of the height 

26^ 



4/r^ at £ = Uj, and secondary, significantly smaller, sharp peaks of the height 4-/[r(r - 1)]^ 
si E - Ui + 0. The situation is illustrated in Figure 13. Therefore, the device constructed 
for fl » 1 works as a branching band-pass filter. Every output j - 2, . . . , r is controlled by a 
dedicated control potential Uj (put on line j + r) which determines the position of the principal 
passband. 
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Figure 13: Characteristics of the branching spectral filter, obtained from the graph in Fig. 12 for « = 2r = 8 and a = 4. 
The control potentials are set to C/2 = 1> C/3 = 0.5, V^ = 0.1. The left figure shows the transmission probabilities 



d(|C/)) 



J^ (E) to the outputs j = 2,3,4, the right figure shows the reflection probability \S\^ {E')\ , the probability of trans- 
mission to the drain and the probabilities of transmission to the controlling lines |S J.| (E)\^ for f = r + 2,. ..,n. 



Identical potentials. If some of the control potentials are of the same nonzero value, or if their 
values are very close, it has effect on the secondary peaks at that energy. Namely, equation (130) 
is changed to 

4 ( he ^ 



■)(U) 



Hm py'(E) ~ -7 ■ — 



E^Ut 



■ht 



for t/f ^ , 



(131) 



where h( can be roughly interpreted as the number of the control potentials that are approxi- 
mately equal to C/f, i.e., h( » # {/ e {2, . . . , r) | a^ a/I - Ui/Ut <K 1 }. Therefore, the height of the 
secondary peak at E - Ue indicates the number of control potentials equal to Uf (or to a certain 
very close value). 

Let us emphasize that the heights of the principal peaks are essentially unafltected. 

The effect of Uj = 0. Setting the control potential Uj on the controlling line j + rtoQ essentially 
results in cancelling the corresponding principal peak. The limit (125) is changed to 



2 1 +/!()+ T 

£^0 Ji I +a^r I +k+ ° 



(132) 
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and, therefore, the principal peak height is changed from the value given by equation (129) to 

limrf\E)^0. (133) 

E—>Uj ■' 

If all the control potentials are zero, then P (E) - for all j - 2, ... , r, and the incoming 
particles are mostly reflected. 

Another operation mode. One may consider the operation mode when the drain line be subjected 
to a nonzero potential Ui. In such a situation a second principal peak at the energy E - Ui 
appears at every output j = 2,...,r. More precisely speaking, if Uj + and satisfies the 
condition a^ ^\l - UqIUj\ » 1, then 

lim !P<.^'(£) » % , (134) 

E^Ui ' r^ 

whereas the relations (127)-(130) are unaffected. Therefore, the device with a nonzero potential 
at the drain line works as a dual-band branching filter with two principal transmission probability 
peaks at the outputs: the first peak is at £ = Uj + 0, where j denotes the output, the second one 
is at £■ = U\. Both peaks are essentially of the same height 4/r^. 

Remark 9.1. The spectral filter designed in Section 6 is obviously a special case of the filter 
studied in this section, corresponding to r = 2. 

10. Practical realization of Fiilop-Tsutsui vertices 

The function of the controllable filter devices designed in the previous sections is based on 
the threshold resonance effects in quantum star graphs with "exotic" couplings in the vertices, 
namely those of the Fiilop-Tsutsui type. It should be emphasized that standard vertex couplings 
(the free and the (^-coupling) would not work that way. It is therefore essential, for the proposed 
designs to be viable, that the required Fiilop-Tsutsui vertices can be physically realized. This 
problem has been addressed in [6] and [7], where it was proved that any Fiilop-Tsutsui coupling 
can be approximately constructed by assembling a few short lines in a web with (^-couplings in 
the vertices and, in some cases, certain vector potentials on the lines. It essentially solves the 
problem, because the (^-couplings themselves have a simple physical interpretation and can be 
well approximated by steep smooth potentials [8]. 

In this section we revisit the results of [6] and [7], proposing a construction which is simpler, 
and, therefore, easier to be experimentally realized. The new construction uses less (^-couplings 
than the scheme from [6], and at the same time, in contrast to the scheme from [7], it requires no 
vector potentials in case that the matrix T from the boundary conditions (20) is real. 

The approximation is constructed in the following way, cf. Figure 14: 

• Take n halflines, each parametrized by x e [0, +cx)), with the endpoints denoted as vy, 

• For certain pairs \j,l] (j,{ e {!,...,«)), join halfline endpoints Vj,ve by a connecting 
line of the length d/yj/, where d is a length parameter and 1 /jj/ is a (positive) length 
coefficient. The pairs {j, {] such that the vertices Vj, V( are connected, as well as the values 
ofyjc, will be specified later If the endpoints Vj, vg are not connected, we may formally set 
y^f - 0, which represents the "infinite length" d/jjj = oo. Note that y^e = yej for all j, (. 
The center of the line connecting Vj and V[ will be denoted by w^c (therefore, Wj/ = wij). 
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Figure 14: An example of the approximation for n = 5. Certain half lines are connected by short lines (in this example 2- 
3, 2-5, . . . ). other are without a direct link (e.g., 1-2, 4-5). The bullets represent 6 potentials, the double line symbolizes 
a vector potential. The connecting lines are generally of different lengths, and each of them carries either a (5-interactions 
(cf. link 1-3), or a vector potential (cf. 1^), or none of them (cf. 1-5, 2-3, 2-5), depending on the parameters of the 
approximated coupling. 



• Put a (5-coupling with the parameter aj(d) at the vertex Vj for each j = 1, ..., n. 

• At certain points Wj/ (to be specified later) place a (5-interaction with a parameter yS;/(J) i^ 
0, at the remaining points Wj/ formally set/3j^e{d) - 0. (Note that Pje(d) = j3ej(d).) 

• On certain connecting lines not carrying a (^-interaction (i.e., those with j3j^e{d) - 0) put a 
constant vector potential supported by the interval \d/(4-yj^(), 3d/(4-yj()\. Its strength will 
be denoted by A(j^{)(d) (if the variable x on the line grows from Vj to V() or A^(j)(d) (if the 
variable x grows from V{ to Vj). Both orientation are allowed, and obviously A^ij){d) — 
-A(j(){d). On the remaining lines formally set Ay (-fid) = 0. 

In the rest of the section we will specify the dependence of jj^, aj{d), /3j^e{d) and A^j^(){d) on 
the matrix T and on the length parameter d. In order to simplify the notation, we will write Uj 
instead of aj{d); the same applies to /3j(, A^j^(y We will also use symbol A^^ standing for the set 
containing indices of all the lines that are joined to the j-th one, i.e. { e Nj <^ jj/ + 0. 

For all y e { 1 , . . . , n), the wave function on the j'-th half line is denoted by the symbol ijjjix), 
X e [0, oo). The wave function on the line connecting the points Vj and vc is denoted by if{j,e), 
X e [0, cZ/yj/]. Obviously, one has to take the line orientation into account. We adopt the 
following convention: If the value corresponds to the endpoint Vj and the value dljj^i to the 
endpoint V(, then the wave function is denoted by fdt). In the case of the opposite direction, the 

wave function is denoted by ifitj). Therefore, ip((j)(x) - ip^^i) I — - x j for all x e 0, — . 

Now we can proceed to determining yy/, aj(d), /3j/{d) and A(j^t){d). At first, let us express 
the effect of the (^-couplings involved. The 5-coupling at the endpoint of the y'-th half line (j - 
1 , . . . , n) implies 



^j{0)^iPud^)^nej)(—\ foralUeA?,., (135a) 

il/p) + Y, ^o;o(0) = aj ■ ilfjiO) , ( 1 35b) 



teNj 
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and the ^-interaction in the point Wj^ means 

"""((41) ='""((41) "'""(4) "•'""''■ "*'* 
'»"((4)J-"«"((4)J=''''-'»"(4)' "'*> 

where (p(x+) and ^(x_) denote the right-sided limit and the left-sided limit, respectively, of (p at 

X. 

Further relations which will help us to find the parameters come from Taylor expansions. 
Consider first the case without the vector potentials (the effect of the potential will be taken into 
account later). We have 



<^) (4;) ^ ^"'^^^^^ ^ 4; ^^'^'^^^ ^ ^^"^'^ ' 
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(138) 



<"((4)J^<'""- 

where the superscript refers to zero vector potential on the short line connecting Vj and V(. If 
the vector potential on a connecting line is nonzero, the relations (137) and (138) are changed 
according to the following lemma. 

Lemma 10.1. Let A eM., s,t e M, m e M\{0). Let i^^, be the solution of the initial value problem 

-^^if,^Aif,, HO)^s, f(0) = f, (139) 

and tf/^, be the solution of the initial value problem 

—- \-ih— - qA{x)\ i//^Ai//, il/{Q) = s, ^'{0) = t , (140) 

2m \ ax 



for a certain q eM. and a function A : R —^ R such that i supp(A), where supp(A) denotes the 
closure of the support of A. Then 

ilf^,(x) = e' i r -^Wdf ■ i/^o ,(x) forallxeR. (141) 

The statement can be proved easily, cf. [15] or [6]. Lemma 10.1 implies that the constant 
vector potential A supported by the interval [Li , L2] shifts the phase of the wave function and of 
its derivative at x > L2 by iqA(L2 - L])lh. Therefore, equations (137) and (138) are transformed 
into 



V(7/) 



— J = e-- ^ J^,.„(0)+_^;.,,(0)J+O(«'2), (142a) 

^o;^)l(2^)j^'"" " ■^o/)(0) + ^(^), (142b) 
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Equations (135), (136), (142) and (143) connect values of the wave functions and their deriva- 
tives at all the vertices Vj, j - I, ..., n. In the next step we eliminate the terms with the "auxiliary" 
functions ^y/). To do so, we at first use equations (136b) and (142b) to obtain 
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Then equation (143 a) gives 
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We substitute for ip(j^()(d/(2jj^()) in equation (145) from equation (142a), and find 
I d \ i^"-^^!! d \ d I d 
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then we apply equation (135a) to rewrite the equation above in the form 

d \ .„s d (^ d 
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hence we obtain 
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and finally we sum the last equation over k, which yields 
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iAf(0)+<9(flf2). (149) 



According to equation (135b), the LHS of equation (149) is equal to -dtl^'XO) + daj\j/j{Q). This 
fact allows us to eliminate the term YIt=\ f'lni^) from equation (149), which leads to 



#;(0) = 
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e~\ 2 + ^fijx 

31 



for all 7 = !,...,«. 

Our objective is to choose aj(d), jj/, /3j^i(d) and AQCf{d) in such a way that the system of 
equations (150) with j = 1, . . . , « is equivalent to the boundary conditions (20) in the limit c/ — > 0. 
It is convenient to adopt the following convention on a shift of the column indices of T: 

Notation 10.2. The lines of the matrix T are indexed from 1 to r, the columns are indexed from 
r + 1 to n. 

Taking this convention into account, we rewrite the boundary conditions (20) as a system of 
n equations, 

n 

•^.(O) + J] f/<A'f(0) = j = 1, . . . , r , (151a) 

r 

0^-J]tijUO) + fj(0) ; = r+ !,...,«. (151b) 



(=\ 



The rows of the system (151) are of two types. We start with the type (151b), corresponding to 
j > r + I. We see immediately that if the parameters aj, jj^, /3j^e and AyQ are chosen such that 






''" " :^ '-ri^^^tj for^<r, (153) 



then (150) tends to 



yj,t^O for€>r+l, (152) 

2 
2 + 

J«, + J]2r,,-2— ^.1, (154) 

r 

= lAi(O) - Yj ^J'Pe^^) for every j^r+\,...,n, (155) 

t=i 

in the limit d/ — > 0, and, therefore, equations (150) for j > r + 1 are equivalent to equations ( 1 5 1 b) . 
It suffices to find a solution of the system of equations (152)-(154). Equation (152) implies that 
there is no connecting line between Vj and v; if ^ > r + 1. If /" < r, we use equation (153); it is 
easy to check that it is satisfied by 

yj,( = \ttj\ , (156a) 

r for tej e M 

^o;0-| j.p,j\^gt,j forr,,ECW ' ^^^^^^ 

^ / forttj>OwtejiR 

Remark 10.3. For every j > r + I, the point Vj is connected to a point V( if and only if ^ < r and, 
moreover, tgj + 0. Furthermore, if v j and vt are connected, then the (5-interaction in w j^c is present 
iff tij < 0, and the vector potential on the short line connecting \j and v^ is present iflF tcj i M. 
Consequently, none of the short lines connected to Vj, j > r + 1, carries both the 5-coupling and 
a vector potential. 
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Finally, we use equation (154) together with (156) to find ay. 
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(157) 



The case y > r + 1 is now solved. Let us proceed to the case j < r. We substitute for i/'^(0) in 
the left-hand side of equation (151a) from equation (150), and in this way we get 
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Then we substitute all the already derived parameters from equations (152), (156) and (157) into 
equation (158). After a manipulation we arrive at 
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Similarly as in the case j > r + 1, we require that the system (159) is equivalent to the sys- 
tem (151a) in the limit c/ — » 0. In this way we obtain a system of (sufficient) conditions on the 
sought parameters Uj, jj( and A(y/) (7 < r), namely 



-ri^^ 2r,;, 



2+ T^Pit u ^, 

2y je"J'^ k=r+l 



+ J] tjkkk = for ^ < r , 



(160) 



1 



1 + 4-^Pj, 1 « 



f=i 



2^/^^.^ ^ 



*=r+l 



2 
5. 



If>*l-i Yj Oi-;i Xi I^^-^I'"^- 



|/:>f|/,i<0) 



/(:=r+l 



(161) 
Hence we derive explicit expressions for aj, jj^t^ f5j^( and A^^) (j < r). From (160) we have 

n 

JiX^ Yj hkkk^{TT)jc, (162) 

therefore, there is a connecting line between v j and v^ if and only if {TT*)j( + 0. Furthermore, 
for (TT*)j( + 0, we use equation (160) to determine j6y/ and A(j/) - 0: 
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If {TT')jc < 0, then 
If {TT*)j( > 0, then 

If {TT*)jt i K, then 



lij,c^Q, A(j,t)^Q. (163a) 



[iji^--{Tr)ji, A(,;,)=0. (163b) 



O Pi 

/3j,e = , Ao;q = -^ ■ - |(rr*),,| ■ arg (-(rr*),,) . (163c) 

Note that on none of the short connecting lines there is both the 5-coupling and a vector potential 
(similarly as in the case y > r + 1, cf. Remark 10.3). 

Finally, we substitute y^e, fij/ and A^j/y found above into equation (161), and obtain the 
expression for Uj, j < r. 



1 



n n r r 

yk=r+\ k=r+\ {k>r\tjt<Q] t=\,ti^i K<''K*jA(7T'),f>0) 



(164) 



Example of the construction 

We conclude this section by a practical demonstration. We use the technique to construct 
approximations of the Fulop-Tsutsui couplings needed for the filtering devices designed in this 
paper. 

Let us consider the coupling used in Section 4, corresponding to T = (a b\. Since r - 1, 
there is no connecting line between the endpoints V2 and V3 due to equation (152), cf. also 
Remark 10.3. If we substitute f^ = a, ti3 - b into equation (156a), we obtain 721 - \a\, 731 = \b\, 
therefore, the endpoints vi and V2 are connected by a short line of the length dl\a\, the points vi 
and V3 by a short line of the length dl\b\, see Figure 15. Since we have assumed a,b eM., we get 
A(3,i) = A(3 2) = from equation (156b); moreover, the device is contructed for a,b > 0, hence 
yS3,i = ^63,2 = due to equation (156c). Therefore, each short line carries neither a ^-interaction, 
nor a vector potential. The parameters of the (5-couplings in the endpoints V2 and V3 follow from 
equation (157): a'2 = (1 - a)/d, Q'3 = (1 - b)/d. Finally, the parameter of the (5-coupling in the 
endpoint vi is given by equation (164): a^ - (a^ + b^ - a - b)/d. The quantity d appearing in 
the formulae is the length parameter. The small size limit (i — > 0, together with the 6 potentials 
parameters properly scaled according to the formulae above, effectively produces the required 
Fiilop-Tsutsui vertex coupling. 

The approximation of the coupling corresponding to T = , used in Section 6, can be 

obtained in the same way. The result is illustrated in Figure 16. 

Finally, let us briefly comment on the case when T = aG for a unitary r x r matrix G; 
the corresponding Fiilop-Tsutsui coupling is essential for constructing the device designed in 
Section 8. It holds: 

• If j, ^ > r + 1, then the endpoints Vj and ve are not connected due to equation (152), 

• If j, € < r, then the endpoints Vj and V( are not connected, which follows from the property 
TT* = a^I and equation (162). 
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^V3 

I •— 

FT. (a b) v^ '^'° V2 



Figure 15: Approximate construction of tlie Fulop-Tsutsui couplings for n = 3 and T = [a by For the special choice 
a = \ corresponding to the filter constructed in Sect. 4 (cf. eq. (36)), the parameters of the S potentials in the points v\ , 
V2 and V3 are q-| = b(b - l)/d, 0-2 = (i.e., there is no (5-interaction in the point V2) and 0-3 = (1 - b)/d. 



^ 




FT: ira) 



Figure 16: Approximate construction of the Fulop-Tsutsui couplings for n = 4 and T = \ . The parameters of the 

r5 potentials in the points V|, i'2, V3, V4 and W2.4 are Q"! = 2a{a- l)/d, 0-2 = 2a(a-2)/d, 03 = (1 - 2a) jd, a^ = {\ -4a)/d 
and/?2,4 = -^ajd, respectively. 

Consequently, the short hnes connect only pairs Vj, vi such that I < r < j. In other words, they 
form a bipartite graph in which {vi, . . . ,vv) and {vv+i, . . . , v„) are independent sets of vertices, see 
Figure 17. In addition, if the matrix G is chosen as a multiple of an Hadamard matrix, then the 
connecting lines are all of identical length dal y/r for d being the length parameter. 



input 
output 




control 



Figure 17: Illustration of the approximate construction of the Fulop-Tsutsui couplings for (i = 2r = 8 and T = aG, 
G e U(r). Couplings of this type are used in Sect. 8. The approximating graph in the center has a bipartite structure. 
The input and the output lines are on the left hand side, the controlling lines are on the right hand side. Some of the 
connecting lines carry a (5-interaction or a vector potential, depending on the phases of tj(. 



11. Scattering matrix of the approximation 

Now we demonstrate the convergence of the approximating scheme proposed in Section 10. 
The proof is performed in terms of the scattering matrices. We calculate the on-shell scattering 
matrix of the approximating graph, and show that it tends to the scattering matrix of the approx- 
imated Fiilop-Tsutsui coupling in the limit d ^> 0. This approach is different from the previous 
papers on this topic [6, 7], where the proof was either omitted, or consisted in performing a full, 
but very lengthy and cumbersome, demonstration of the norm-resolvent convergence. In addi- 
tion, the convergence of the scattering matrices is proved here in a general setting when the lines 
coupled in the vertex carry potentials Vi, . . . , y„, which is another novel point of this paper. 
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In what follows, we denote the on-shell scattering matrix of the approximating graph by 

SdiE^Vu ■ ■ ■ ,yn) ■ (165) 

The matrix Sd(E; Vi, . . . , V„) will be calculated in the following three steps: 

1. We derive relations between the on-shell boundary values i//j{0) and derivatives t//':(Q) by 
analyzing the physical properties of the internal web. 

2. We substitute the final-state wave function components *P (x), expressed in terms of 
Sd(E; Vi, . . . , V„) (cf eq. (10)), into the system of equations obtained in the first step. 

3. We solve the system, in this way we find Sd(E; Vi, . . . , V„). 

Step 1. Let tl/j and ipje (j - I, . . . ,n, { e Nj) be the components of a wave function on the 
approximated graph, corresponding to a given energy E. They satisfy the boundary conditions 
in the vertices; our goal is to derive relations among the boundary values lA/O), i/''(0). The 
5-coupling at the point Vj implies 



2 <f'j,iO) + ifr'jiO) = ajil^jiO) for all j = 1 , 



. ,n , 



(166) 



heN. 



where <Pjh(x) is the wave function component on the line connecting Vj and v/,, with the conven- 
tion X = at Vj and x = d/yj^h at v;,. 

We want to express f'/jiO) in terms of iA>(0) and t//i,(Oy, the result will be then substituted 
into equation (166). Let us assume at first that the short line connecting Vj and v/, carries a 6- 
interaction with parameter /3j^h in the center The corresponding wave function component is 
given by 

C+e'*-' -H C^e-'*-' for all 



iPjhix) = 



j^+^ikx ^ £,-g 



-ikx 



xe 



for all X e 



k/ 1 

' 2ri./, 


' 


d a 

T-y,.h ' y, 


h 



(167) 



where C^,C ,D^,D are certain constants and 



2mE 



(168) 



Everywhere in this section, the symbol k stands exclusively for the quantity given by equa- 
tion (168). From equation (167), we have 



ifrjiO) ^ipjhiO) ^C^ + C-, 
d 



i/^h(0) ^fjh 



7j,h 



= D+e ''.'.» + D e ''/.» 



(169a) 
(169b) 



Furthermore, the (^-interaction in the center of the line means 

and 



d d 

fjhii;: )-) = ^j7i((:^— )+) 



2rj,h 



2yj,h 



"f'jhii^^—)^) - "P'jhii^^Z-:^-^ ^P},hVjh{w—) , 



2yj,h 



2yj,h 



2yj,h 



hence 



C+e'*'''/'' +C'e '*'''/.'■ = D+e'* '''-'■■'■ + De '*'''J.'" 



i/tC^e -'''■'' -ifeC^e ''J'- - i/tD+e ''''■'' +ikD e 
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-h^(c'^'' 



(169c) 

^ -HC^e"'*'^] . (169d) 



By solving the system (169), we find C^ , C and D^,D in terms of i//j{0) and i//h(0)- Then we 
substitute C^, C into the equation 

ip'j,(0)^ik{C^-C-), (170) 

which follows from (167), and in this way we find 

f'jhiO) = ^^ ^— i-ifrjiO)Lsk— + ^ sinfe— ) + MO)] . (171) 

yjj, k 2yjj, 

Recall that the formula (171) has been derived on the assumption that the short line connecting 
Vj and Vfi carries a ^-interaction with the parameter jSj/,. The solution for the case when the line 
carries a vector potential A(y/,) can be obtained using Lemma 10. 1 . Both solutions can be merged 
into the general formula 

^./.(O) - ;, 1^ ^-r (-"^/O) f'^°s k—+^ sink—] + e"'4r^<^V,(0)) (172) 

^ sinfe^ + ^sin^/t^ I I 7j,h 2fe yjjj J 

7j.h k Ijj^k 

which is applicable to any of the three possible situations (a (5-interaction, a vector potential, 
none of them), it suffices to set/3j^h = or Ay^h) = if necessary. Now we sum equations (172) 
over h and apply (166), in this way we get n equations relating i//j(0) and (^',(0): 

-^ E 7^. T^ ■'A/0)+fe y J- ---^,(0) + (A;(0) = ajifrj(O) , 



^pw, sinfc-^ + ^ sin ^T^- i^ sinfc-^ + ^ sin fc^ — 

"'="1 yj.h k 2yj.h '"=™./ yj.h k 2y^^ 



;=l,...,n. (173) 



Step 2. If a particle with energy E comes from the ^-th line into the internal approximating 
web, the final-state wave function naturally obeys relations (173). Let ^F. (x) denote the com- 

ponents of the final-state wave function on the half lines and *P (x), (*P ) (0) be the boundary 
values and derivatives. We express *Pj(0) and *I".(0) in terms of the entries of Sd{E; Vi, . . . , V„) 
using the definition of the scattering matrix (10): 



i/W/ 



-j^(l + [Sd]a) for j^i. 



n^(0).j ^'_L^[^,]^, fo,,-^,, (174) 



(^^^^(O)^^^^"^'^"^ f"':^' (1^5) 



] 



where kj - ^2m{E - Vj)/h?- for all j = 1 , . . . , n. 

If we multiply equation (173) by d, then for all y = 1, . . . , n replace iffj(0) and lA (0) by ^ . (0) 
and (¥.^^)'(0), respectively, and finally, substitute for ¥^\Q) and (¥^V(0) from equations (174) 
and (175), we obtain a system of n^ equations for [Sd]je, j,( - 1, . . . ,n. They take different 
forms for j i^ { and j - i: 
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Case j + 



heN: 



cosk— + 



Pi.li 



sinfe- 



7/,/. ■ 24 y,;A 

^ aid + xkjd 



kdY^ 



[Sd]jt+ 



e -''.'■'■ 



teWf 



sinfe^ +^sin2/t^ Vfe^ 



[^rf: 



/i^ 



2r;,/. 



'2777 ''<-'■■'■' 



-M 



Ji.f k T-y,i.c 



■XN,(i), (176a) 



wheiexN (0 - 1 for ^ e A^^ andxN (i) - otherwise. 
• Case j = { 



cosfc^+^sinfe^ 
-kd } , ^ — J ;; 7 a{d + iked 



i^.sinfe^+^sin^^^ 



heN, 



'iye 

■kdYj 

heN, 



yke 



'l^^i'M 



sink^ + ^-^sin^k:^ ^Ih 

yt.h k 2y,,i 



[Sd] 



he 



._ cos A:— 



^sinfe^ 



AeAff 



sinfe^ + q^sin^fe:^ 



+ CK^C/ + ifef c/ 



ytj, 



T-yt,h 



1 



(176b) 



Step 3. The sought scattering matrix Sd(E; Vi, . . . , V„) can be obtained by solving the system 
of equations (176). We observe that the system (176) can be rewritten in the matrix form 



(ZD ' + idD)Sd = -ZD ' + idD, 
where D = diag( yfici, . . . , yfkn) is the matrix introduced by equation (15) and 



(177) 



'-jt 



kd- 



■XN,{{) 



for j ^{, 



yj.c 



2r// 



-kdY, 



heNi 






Tf,;, ' 27f ft 



(178) 



Therefore, 

SdiE; Vi, . . . , y„) = (ZD"' + idDy\-ZD-^ + idD) = -/ + 2idD(z + idD^^ ^- (1^9) 

Remark 11.1. Note that the matrices Z and D occuring in equation (179) have very different 
meanings. The matrix Z represents primarily the internal properties of the approximating graph, 
and depends on the approximated boundary conditions, on d and on E. The matrix D represents 
only the external setting, and depends on Vi , . . . , y„ and on E. 

It suffices to determine Z. This will be done by substituting the values of the coupling param- 
eters Uj, Pjx, the length factors y^t and the vector potential strengths A(^/) into equation (178). 

38 



All these parameters are given by equations (156), (157), (162), (163), (164). We shall keep 
in mind that the expressions are different when j, { belong to { 1, . . . , r) and when j, { belong to 
{r + 1, . . . , «), and moreover, Zjc is defined in different ways for j i^ { and for j = {. Therefore, 
the formula for Zjc is divided into several expressions: 

• Case j,{ < r and j ^ ^ 



Zje- 



-We'"g<^^'>'f(sinj(54%;i) 



if (TT*)j( < or (TT*)je ( R , 



-Me'"''s(^^*'^Msin 



kd 



\(TT')jc\ 



;I(^?")/H.sjjj2 kd 







Case j - { < r 



Zte - - 



Z 



kd cotg 



kd 



kd 



2\(TT 



^r if(rr),,>o, 
if(rr*),, = 0. 



(180a) 



{h<rM{\(TT')i,e<OvtMi« 



liTT^M 



Z 



cos 



kd 



kd- 



\{TT'),f\ 



a\(JTMc\„ kd 
■ ^ kd ^"' 



\(TT'h,\ 



kd 



{h<rMl\(TT')u<Q) 



\(TT'%, 



> \(TT')i,f\ -2 kd 
' kd *"' 2\(TT')u\ 



( n 



h=r+\ (A>r|ta<0) {h<rMe] {h<r\hi=(A(TT')tk>0} 



\h=r+l 



(180b) 



Case j < r < { 



'-jt 



kde''^'^ 



''i-"R^r 



if tj( > or tji i : 



Me'^''g''Msin 



'i' (s 



kd_ _ oM 





sm-^l ifO^<0, 



(180c) 



2kxl / 



iff;^ = 0. 



Case { < r < j 

Case j,{ > r and j ^ ^ 
Case j - { > r 



where Zcj is given by equation (180c) . 



Z,-, = 0. 



(180d) 
(180e) 



Zff — - 



^ M ^ cos 1^ - 4!|^ sin ^, 

> fed cotg ■ - > kd — T—. ;; — 

Zj ^ \t,.f\ Z-i sin M -sl^cin^ -^ 



1 - z \^h{\ + Z ^^''^ 

{h<i]tht<0] 



\ h=\ 



(180f) 



Equations (180) above fully determine the matrix Z. Therefore, for a given T, V\, . . . ,V„, E 
and d, we can easily find SdiE; Vi, . . . , V„) using the formula (179) together with equations (11), 
(15) and (168). 
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If d — » 0, i.e., kd «: 1, we can expand each of the expressions (180). The calculation leads to 

z^i^'ir _^ji_,) j+o((M)2). (181) 

When the matrix (181) is substituted into equation (179), one obtains the following approxima- 
tive formula for Sd{E; Vi, . . . , ¥„): 

S,{E- Vi, . . . , y„) = -/ + 2 ( gS^w . j (efi) + TQl;r*\' ( gd) TQ^^) ) + 

^H G(2)f/-rw-^rey )( ef2)r*w-'e(,) {i-Q',,,rw-^T)Q,,, ) + 0(d'), (182) 

where W - Q}^. + TQ%T* and 2(i), 2(2) are the diagonal matrices introduced in Proposition 3.1. 
In equation (182), the scattering matrix of the approximating graph is found up to a small error. 
The error tends to zero for d -^ 0, i.e., 

lim>S^(£; Vi, . . . , y„) = -/ + 2 ( ^^^^, j {q^^ + TQl^rf' ( 6(1) rgfz) ) . (183) 

Now we compare the limit (183) with the scattering matrix of the approximated star graph. Since 
the vertex of the approximated graph carries a Fulop-Tsutsui coupling, described by the bound- 
ary conditions (20), its scattering matrix S(E; Vi,. . . , ¥„) is given by formula (25). Comparing 
equations (183) and (25), we find 

limSdiE; Vu-.., V„) = S(E; Vi, . . . , y„) . (184) 

To sum up, the scattering matrix of the approximating system converges to the scattering matrix 
of the requested Fiilop-Tsutsui vertex coupling in the limit c/ — > 0. 

12. Conclusions 

The resonance phenomenon found in this paper allows to control the transmission of a quan- 
tum particle along a line by external potentials put on other lines connected to the channel at a 
vertex. We have demonstrated that for certain well-chosen Fiilop-Tsutsui vertices, the resonance 
is strong enough to enable a design of quantum spectral filters with fine characteristics. The scat- 
tering characteristics are very different depending on the matrix T determining the Fiilop-Tsutsui 
coupling in the graph vertex. Therefore, it is likely that many more controllable quantum devices 
based on this concept can be obtained for other choices of T. 

Some of the designed filters have rather small transmission probabilities inside the passbands. 
It would be useful to find how to increase the probabilities, and in this way to improve the filter 
characteristics. 

Macroscopic controllability can be achieved also by an application of a magnetic field, for 
instance, by tuning a magnetic flux through rings carried by the controlling lines, or through 
ring-shaped controlling lines. It turns out that the use of suitably chosen Fiilop-Tsutsui couplings 
together with the application of magnetic field lead to transmission characteristics with adjustable 
high peaks. We plan to address this problem in a subsequent article. 
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It seems that the idea can be extended to the control of the spin flow, because the mathemat- 
ical description of a spin one-half particle propagating in a quantum star graph with n lines is 
equivalent to the description of a spinless particle in a quantum star graph with 2n lines. Explor- 
ing the spin-filtering mechanism could be of interest. 

Finally, there are known parallels between the particle propagation in quantum graphs and the 
wave propagation in waveguides. Consequently, we are convinced that the results of this paper 
can be adapted to the control of a transmission of macroscopic electromagnetic microwaves in a 
thin waveguide. 
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